LEI^ 


ON  THE  FIRST  TIME  A SEPARATELY  MAINTAINED 
PARALLEL  SYSTEM  HAS  BEEN  DOWN  FOR  A FIXED  TIME 


by 

SHELDON  M.  ROSS 
and 

JACK  SCHECHTMAN 


» LU 


OPERATIONS 
RESEARCH 
CENTER  , 


ORC  78-5 
APRIL  1978 

r # — _ X 


DOC 


i l ) 

SEP  2T  1978  l!i 


r 


2flED 

D 


I9T8  ’ 

nnai 


Approved  for  public  release; 
Distribution  Unlimited 


UNIVERSITY  OF  CALIFORNIA  - BERKELEY 


tNIA  -BERKELEY  _ . 

78  08  16  046 


ON  THE  FIRST  TIME  A SEPARATELY  MAINTAINED  PARALLEL 
SYSTEM  HAS  BEEN  DOWN  FOR  A FIXED  TIME 


by 


Sheldon  M.  Ross 

Department  of  Industrial  Engineering 
and  Operations  Research 
University  of  California,  Berkeley 

and 

Jack  Schechtman 

Institute  of  Pure  and  Applied  Mathematics 
Rio  de  Janeiro,  Brazil 


D D C 


aEL'  •'] 

SEP  W 1978  : 


EEGDI/lslllj 

D 


APRIL  1978  ORC  78-5 

This  research  has  been  partially  supported  by  the  Air  Force  Office  of 
Scientific  Research  (AFSC),  USAF,  under  Grant  AFOSR-77-3213  and  the 
Office  of  Naval  Research  under  Contract  N00014-77-C-0299  with  the 
University  of  California.  Reproduction  in  whole  or  in  part  is  per- 
mitted for  any  purpose  of  the  United  States  Government. 

~DlSTRIBUnON  <-  | 

Approved  for  puU;.:  . 

Di»tri  bution 


Unclassified 

SECURITY  CLASSIFICATION  of  THIS  FACE  (Whin  Date  Entered) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


2.  GOVT  ACCESSION  NO 


1.  RE 


ENT'S  CATALOG  NUMBER 


£ 


RE  OF  REPORT  S-RERIOO  COVERED 


ON  THE  FIRST  JIME  A SEPARATELY MAINTAINED 
PARALLEl  SYSTEM  HAS  "BEEN  DOWN  FOR  A FIXED  TIME 


/ Research  Xeport 


! 


t.  RERFORMING  PRO.  REPORT  nTmTeR 


7.  AUTHOR^ 

S - * ! 

Sheldon  My/Ross  and  Jacky^chechtman  , 


• CONTRACT  OR  GRANT  NUMBERff) 

i/AFOSR-77-3213  • 


16.  rrogram  element,  rroject,  TASK 
AREA  * BBRK.  UNIT  NUMBERS 

jy 

'0>J  23J04//A5 


» rerformino  organization  name  ano  ADDRESS 
Operations  Research  Center 
University  of  California 
Berkeley,  California  94720 


W 


<1.  CONTROLLING  office  name  ano  aooress 

United  States  Air  Force 
Air  Force  Office  of  Scientific  Research 
Bolling  AFB,  D.C.  20332 


*1  flBRBRT  PATE 

Jtortt  »78/ 


IS.  NUMBER  OF  RAOES 

15 


1 4.  MONIXORING  AGENCY  name  A AOORESSCI/  dlllerent  Inm  Controlling  Olllei) 


IS.  SECURITY  CLASS,  f of  Uilt  ripen) ■ 

Unclassified 


Ua.  OBCLASSI  FI  CATION/ DOWN  G R ADl  NO 
SCHEDULE 


l«.  DISTRIBUTION  STATEMENT  ( o I Ihii  Report) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (o I Ihi  ibitreet  intend  In  Block  JO,  II  dlllerent  Item  Ripen) 


II.  SURRLEMENTARY  NOTES 


Also  supported  by  the  Office  of  Naval  Research  under" CStitrdct 

mvu-n- c-fuy  ( 


IS.  KEY  WOROS  (Continue  Si  reverie  ildi  II  nicimrr  wtd  Identity  by  block  number) 


Parallel  System 
New  Better  Than  Used 
Down  for  Fixed  Time 
First  Passage 


20.  AES!  RACT  (Continue  on  reverie  ildi  II  memory  end  Identity  Ay  Slock  number) 


(SEE  ABSTRACT) 


DD 


FORM 
I JAN  71 


1473  EDITION  OF  I NOV  ft  IS  OBSOLETE 

S/N  0102-  LF-  014-  6601 


Unclassified 


J//0  r/C0 


SECURITY  CLASSIFICATION  OF  THIS  RAGE  (Whin  Dele  tnh 


ABSTRACT 


Consider  a system  consisting  of  n separately  maintained  independent 
components  where  the  components  alternate  between  intervals  in  which 
they  are  "up"  and  in  which  they  are  "down."  When  the  iC^  component 
goes  up  [down]  then,  independent  of  the  past,  it  remains  up  [down] 
for  a random  length  of  time  having  distribution  F^[G^]  and  then  goes 
down  [up].  We  say  that  component  i is  failed  at  time  t if  it  has 
been  "down"  at  all  time  points  s e [t  - A,t]  ; otherwise  it  is  said 
to  be  working.  Thus  a component  is  failed  if  it  is  down  and  has  been 
down  for  the  previous  A time  units.  Assuming  that  all  components 
initially  start  "up"  let  T denote  the  first  time  they  are  all  failed, 
at  which  point  we  say  the  system  is  failed.  We  obtain  the  moment 
generating  function  of  T when  n ■ 1 , for  general  F and  G , thus 
generalizing  previous  results  which  assumed  that  at  least  one  of  these 
distributions  be  exponential.  In  addition  we  present  a condition  under 
which  T is  an  NBU  (new  better  than  used)  random  variable.  Finally 
we  assume  that  all  the  up  and  down  distributions  F^  , , i * 1,  . . . , n , 

are  exponential  and  we  obtain  an  exact  expression  for  E(T)  for  general 
n ; in  addition  we  obtain  bounds  for  all  higher  moments  of  T by  showing 


that  T is  NBU. 


ON  THE  FIRST  TIME  A SEPARATELY  MAINTAINED  PARALLEL 
SYSTEM  HAS  BEEN  DOWN  FOR  A FIXED  TIME 
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by 

Sheldon  M.  Ross  and  Jack  Schechtman 

0,  INTRODUCTION  AND  SUMMARY 

\ 

-/'in  considering  a system  that  works  for  a random  time  and  when 
failed  is  fixed  in  a length  of  time  that  is  also  random  an  important 
question  is  the  study  of  the  first  time  the  system  is  not  working  for 
an  interval  of  time  longer  than  some  prespecified  value.  For  instance 
in  a nuclear  reactor,  when  the  safety  system  is  out  for  some  critical 
time  it  is  necessary  to  shut  down  the  complete  system  with  all  the 
problems  this  entails.  In  the  food  industry  where  food  must  in  general 
be  kept  at  a certain  temperature,  an  important  question  when  the 
refrigeration  system  goes  down  is  how  long  this  situation  can  be  main- 
tained before  the  food  becomes  spoiled. 

In  "this  paper' we 'consider4  a system  consisting  of  n separately 
maintained  independent  components  where  the  components  alternate  between 
intervals  in  which  they  are  "up"  and  in  which  they  are  "down."-  When 
the  iC^  component  goes  up  [down]  then,  independent  of  the  past,  it 
remains  up  (down]  for  a random  length  of  time  having  distribution 
F^[G^]  and  then  goes  down  [up].  We  say  that  component  i is  failed 
at  time  t if  it  has  been  "down"  at  all  time  points  s e [t  - A,t]  ; 
otherwise  it  is  said  to  be  working.  Thus  a component  is  failed  if  it 
is  down  and  has  been  down  for  the  previous  A time  units.  Assuming 
that  all  components  initially  start  "up"  let  T denote  the  first 
time  they  are  all  failed,  at  which  point  we  say  the  system  is  failed. 
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In  Section  1 we  obtain  the  moment  generating  function  of  T when 
n ■ 1 , for  general  F and  G , thus  generalizing  results  in  [2]  and 
[3]  which  assumed  that  at  least  one  of  these  distributions  be  exponential. 

In  Section  2 we  present  a condition  under  which  T is  an  NBU  (new  better 
than  used)  random  variable.  In  Section  3 we  assume  that  all  the  up 
and  down  distributions  , i * 1,  . . . , n , are  exponential  and 

we  obtain  an  exact  expression  for  E(T)  for  general  n ; in  addition 
we  obtain  bounds  for  all  higher  moments  of  T by  showing  that  T is  NBU. 
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1.  The  Case  n ■ 1 


Let  us  denote  by  N the  number  of  "up"  intervals  that  occur  before 
the  component  fails.  Then  given  N ■ k , we  can  represent  T by 


T - X.  + •••  + X.  + + •••  + Yf  . + A 

1 k 1 k-1 


th  A 

where  denotes  the  length  of  the  i up  cycle  and  the  length  of 

the  i^  down  cycle.  All  the  random  variables  in  the  representation  (1) 
are  independent  with  the  X^  having  distribution  F and  the  Y^  having 


distribution 


'K  - - 


P{Y  < x Y < A} 


< Al  * l 


x < A 


x > A 


where  F is  the  distribution  of  an  up  cycle  and  G that  of  a down  cycle. 


P{N  - k}  - G(A)(G(A))  , k - 1,  ... 


where  G ■ 1 - G , we  obtain  the  moment  generating  function  of  T by 
conditioning  on  N as  follows. 


E[e8T] 


- E[E[e“ 1 j Nl  ] 


‘ e[*SA<*x<s))N(^a<s))N’1 


eS%x(s)G(A)  l ^x(s)4>  A(s)G(A) 


r 


eS%x(s)G(A) 

1 - G(AHx(s)$  A(s) 
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ii 


\ 


where 


t 


^(s)  * E[esX] 


* A(S) 


go 

* J esxdF(x) 

0 

E[eSY  ' Y < A] 


f esxdG(x) 

J G(A) 


For  the  special  case  in  which  X is  exponential  with  mean  1/A 


and  Y is  exponential  with  mean  1/u  we  have 


[>T]  . Hu  - s)e~(|J~s 

^ ^ - (A  + p)s  + Ape  e 


a result  previously  obtained  in  [2]  and  [3]. 

All  of  the  moments  can  now  be  obtained  by  successive  differentiation 
of  (2),  or  by  a direct  conditioning  argument.  For  instance  we  obtain 


(3) 


E[T]  - E[E[T  | N]  ] 

- E[NE[X]  + (N  - 1)E[Y  | Y < A]  + A] 


. mi 

G(A) 


n 

/ 


xdG(x) 


G(A) 


+ A . 


By  viewing  the  working-failed  system  as  an  alternating  renewal 
process  it  follows  that  the  long  run  proportion  of  time  the  component 


i 


is  failed  is 
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2.  WHEN  IS  T NBU , n « 1 

The  nonnegative  random  variable  W is  said  to  be  new  better 
than  used  (written  NBU)  if 

P{W  > s + t j W > s}  <_  P{W  > t}  Vs  , t >_  0 . 

If  we  think  of  W as  representing  the  life  of  some  object  then  W 

NBU  means  that  the  additional  remaining  life  of  any  s year  old  (i.e., 

used)  item  is  stochastically  smaller  than  that  of  a new  item,  for  all  s . 

If  W is  NBU  and  has  distribution  function  H then  we  also  say 
that  H is  NBU. 

Proposition  1: 

If  X , the  length  of  an  up  time,  is  NBU  then  so  is  T . 

Proof : 

Suppose  failure  has  not  yet  occurred  by  time  s . Now  there  are 
2 possibilities: 

Case  1: 

At  time  s the  component  is  up  and  has  been  up  for  a time  t . 

In  this  case  the  remaining  time  to  failure  has  the  distribution  of  the 
convolution  of  Ft  and  H , where  F is  the  distribution  of  remaining 
up  time  for  a component  that  has  been  up  for  a time  t and  H is  the 
distribution  of  time  to  failure  starting  with  the  component  initially 
down.  But  since  is  stochastically  smaller  than  F (definition 

of  X being  NBU)  this  distribution  is  stochastically  smaller  than  the 
convolution  of  F and  H , which  is  the  distribution  of  T . 
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Case  2: 

At  time  s the  component  is  down  and  has  been  down  for  a time  t 
(necessarily,  t * A) . In  this  case  the  remaining  time  to  failure  has 
some  distribution  call  it  D . However  the  distribution  of  T can  be 
written  as  the  convolution  of  D and  the  distribution  of  the  first 
time  that  the  component  has  been  down  for  t consecutive  time  units. 
This  latter  convolution  distribution  is  clearly  stochastically  larger 
than  D . 

Thus  in  all  cases  the  distribution  of  T is  stochastically  larger 
than  the  distribution  of  remaining  time  until  failure.  Hence  T is 
NBU . | | 


3.  EXPONENTIAL  LIFETIMES,  GENERAL  n 


In  this  section  we  suppose  there  are  n components  and  the 
distribution  of  up  [down]  time  for  the  i^  component  is  exponential  with 
rate  >.  [y^]  , i - 1,  ....  n . We  start  ay  deriving  E { T ] , the  expected 
time  until  the  system  fails,  that  is  until  all  components  are  failed, 
starting  with  all  components  up. 

We  can  write  T as  the  sum  of  independent  random  variables  as 
follows 


(4) 


T - T 


A-0 


+ Z 


where  T,  _ denotes  the  first  time  that  all  components  are  down  (it  is 
A-U 

thus  equal  to  T in  the  special  case  A-0)  and  Z the  extra  (or 
additional)  time  from  TAmQ  until  all  components  are  failed.  Now  Brown 
in  [1]  has  computed  E ] and  showed  it  to  equal 


■ l l 


k Ui 


n r-L  - (-1)1 
j-1  \ 


k-1  i < i <• • • <i 
12  k 


l l\  * ui  ) 

J-1  \ i if 


Thus  it  remains  to  compute  E [ Z ] . Let  M denote  an  exponential  random 


variable  with  rate  y h \ y.  . Then  by  conditioning  on  whether  or  not 

1 

all  components  remain  down  in  the  A time  units  following  time  we 

obtain 


E [Z ] - Ae~uA  + (1  - e“uA)[E[M  | M A]  + E[DJ  + E[ZJ] 


9 


where  D Is  the  time  until  all  components  are  down  given  that  they 


were  all  down  and  one  has  Just  gone  up.  Thus,  from  the  above,  we 
obtain 


“MX  , 

pxe  dx 


E(Z]  - A + (eUA  - 1)  — 


+ E[D] 


1 - e 


However,  Ross  in  [5]  has  shown  that 


E[D]  - 


1 “ n TT~ 

J 

n n A 

l n ^r~hr 

1-1  J-i  ui  J 


and  thus  the  expression  for  E[l]  follows  from  (A),  (5)  and  (6). 

The  next  proposition  partly  characterizes  the  distribution  of 
T and  will  enable  us  to  obtain  bounds  on  all  higher  moments  of  T 


Proposition  2: 


T is  NBU. 


Proof : 


Suppose  that  all  components  have  never  been  simultaneously  failed 
by  time  s . There  are  2 cases. 


Case  1: 


At  time  s all  components  are  down,  the  one  that  has  been  down 
for  the  shortest  time  having  been  down  for  a time  t (where  necessarily 
t < A) . Since  T can  be  expressed  as  T.  (the  first  time  all 

r\  u 


1 
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components  have  been  down  for  the  past  t time  units)  pius  a random 
variable  having  the  same  distribution  as  the  remaining  time  to  failure 
of  the  system,  it  follows  that  T is  stochastically  larger  than  the 
remaining  time  to  system  failure  in  this  case. 

Case  2: 

Not  all  components  are  down  at  time  s . In  this  case  the  remaining 
time  to  system  failure  can  be  written  as  the  time  until  all  components 
are  down  plus  an  independent  random  variable  having  the  same  distribu- 
tion as  Z in  the  representation  (4).  Now  Ross  in  [6]  has  shown  that 
the  time  until  all  components  are  down  is  stochastically  larger  starting 
with  all  being  initially  up  than  starting  in  any  other  position.  Hence, 
from  the  representation  (4),  it  follows  that  the  remaining  time  to 
system  failure  at  time  s is  stochastically  smaller  than  T . 

Hence,  in  all  cases  T is  stochastically  larger  than  the  remaining 
time  to  system  failure;  thus  proving  the  result. || 

The  above  result  is  particularly  useful  as  it  enables  us  to  obtain 
bounds  on  E(f(T)]  whenever  f is  an  increasing  convex  function,  by 
use  of  the  following  special  case  of  Theorem  4.6  of  Marshall  and  Proschan 

fM  • 

Proposition  3: 

If  X is  NBU  with  mean  1 / A , then 

no 

J f(x)Ae  *XdA 
0 


E[ f (X) ] 


< 


for  all  Increasing  convex  functions  f . 

In  words  Proposition  3 says  that  if  X is  NBU  then 
E I f (X ) ] F[f(M)]  for  all  increasing  convex  f , where  M Is  an 

exponential  random  variable  having  the  same  mean  as  X . 

Corollary  1: 


Var(T)  < ( E ( T ] ) 2 . 


Proof : 

Follows  immediately  from  Propositions  2 and  3 by  use  of  the 
2 

function  f(x)  * x 


12 


REFERENCES 


fl]  Brown,  M.,  "First  Passage  Distribution  for  Exponential  Repairable 
Parallel  Systems, " Berkeley  Symposium  on  Reliability  and 
Fault  Trees,  SIAM  publications,  (1975). 

[2]  Von  Ellenrieder,  A.  and  A.  Levine,  "The  Probability  of  an  Excessive 

Non-Functioning  Interval,"  Operations  Research,  Vol.  14, 

No.  5,  (1966).  

[3]  Von  Ellenrieder,  A.  and  J.  Schechtman,  "Sobre  El  Problema  Del 

Intervalo  Critico  de  No  Functionamento,"  Atas  do  I Simposio 
Brasileiro  de  Pesquisa  Operacional  e suas  Aplicaijoes,  Vol.  II, 
pp.  534-541,  Brasil,  (1968),  (in  Portuguese). 

[4]  Marshall,  A.  W.  and  F.  Proschan,  "Classes  of  Distributions 

Applicable  in  Replacement,  with  Renewal  Theory  Implications," 
Proceedings  of  the  6th  Berkeley  Symposium  in  Mathematical 
Statistics  and  Probability,  University  of  California  Press, 
Vol.  1,  pp.  395-415,  (1972). 

[5]  Ross,  S.,  "Asymptotic  Properties  of  Multicomponent  Systems," 

Berkeley  Symposium  on  Reliability  and  Fault  Trees,  SIAM 
publication,  (1975). 

[6]  Ross,  S.,  "On  Time  to  First  Failure  in  Multicomponent  Exponential 

Reliability  Systems,"  Journal  of  Stochastic  Processes  and  its 
Applications , Vol.  4,  pp.  167-173,  (1976). 


- ( 1—  I mmm  — 

- . — : -At,.  . . d. 


